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supercondnctors 
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The gapless nodal fermions exhibit non-Fermi liquid behaviors at the nematic quantum critical 
point that is supposed to exist in some d-wave cuprate superconductors. This non-Fermi liquid 
state may be turned into a disorder-dominated diffusive metal if the fermions also couple to a 
disordered potential that generates a relevant perturbation in the sense of renormalization group 
theory. It is therefore necessary to examine whether a specific disorder is relevant or not. We study 
the interplay between critical nematic fluctuation and random chemical potential by performing 
renormalization group analysis. The parameter that characterizes the strength of random chemical 
potential is marginal at the one-loop level, but becomes marginally relevant after including the 
two-loop corrections. Thus even weak random chemical potential leads to diffusive motion of nodal 
fermions and the significantly critical behaviors of physical implications, since the strength flows 
eventually to large values at low energies. 

PACS numbers: 73.43.Nq, 74.72.-h, 74.25.Dw 


I. INTRODUCTION 

The low-lying elementary excitations of d-wave cuprate 
superconductors are known to be massless nodal fermions 
that have a linear dispersion and fulfill the relativistic 
Dirac equation . These fermions govern many of the un¬ 
usual low temperature properties of the superconducting 
state. In the past twenty years, there have been exten¬ 
sive experimental signatures [Tj supporting the fact that 
the nodal fermions are nearly non-interacting and have 
a rather long lifetime. However, this feature can be sig¬ 
nificantly changed if the fermions interact with certain 
critical bosonic mode. Recently, various measurements 
observed an anisotropy in the physical properties of some 
cuprate superconductors Such an anisotropy is 

usually attributed to the emergence of a novel electronic 
nematic order [m, which spontaneously breaks C4 
symmetry of the system down to C2 symmetry. Both 
experimental 0-0 and theoretical studies sug¬ 

gest that a zero temperature nematic quantum critical 
point is therefore expected to exist somewhere in the su¬ 
perconducting dome. In the vicinity of this point, the ne¬ 
matic order parameter fluctuates quantum-mechanically 
around its vanishing mean value. The critical nematic 
fluctuation couples strongly to the nodal fermions, which 
gives rise to severe fermion damping [Tsl - I^ and other 
striking properties (23-1^. 

The quantum critical phenomena associated with the 
critical nematic fluctuation become more interesting, and 
meanwhile more complicated, when there are certain 
amount of disorders. Disorders play an essential role 
in modern condensed matter physics [3l|, |32j, and can 
result in a plenty of prominent phenomena, such as An¬ 
derson localization and metal-insulator transition. The 
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past two decades have witnessed intense research activ¬ 
ities devoted to the study of two-dimensional massless 
Dirac fermions moving in a disordered potential. Realis¬ 
tic systems that contain these fermions include the afore¬ 
mentioned d-wave cuprate superconducto rs |T^ . [s^ [s^ , 
graphene [33 - 1^ , quantum Hall effect [40l - l42j| . and topo¬ 
logical insulators [^. In systems composed of Dirac 
fermions, there can be three sorts of disorders [d^. Ac¬ 
cording to the coupling with nodal fermions, the disor¬ 
ders might be random chemical potential, random mass, 
or random gauge potential [d^ . Theoretical analysis has 
demonstrated that these disorders can produce different 
behaviors of Dirac fermions [TsI . [33 - ld3| . 

Back to the d-wave cuprate superconductors, an im¬ 
portant question is whether these disorders drive an in¬ 
stability of the nematic quantum critical point. It is also 
interesting to examine the influence of disorders on the 
low-energy properties of nodal fermions. Recently, we 
have studied this problem by means of renormalization 
group (RG) techniques at the one-loop level [H, [d^ . In 
particular, we have introduced a parameter C to charac¬ 
terize the effective strength of disorder, and calculated 
the RG flows of C after taking into account the inter¬ 
play between fermion-nematic interaction and fermion- 
disorder interaction. It was found that C, vanishes in 
the low-energy region in the presence of random mass 
and random gauge potential |d6l |. However, C does not 
flow at all in the case of random chemical potential, 
i.e., dC,ldl = 0 with I being a running length scale (d3| . 
Therefore, the parameter ^ for random chemical poten¬ 
tial is considered as marginal. In the other two cases, 
dC/dl < 0, so C vanishes in the low-energy regime and is 
hence irrelevant. If the perturbative expansion is reliable, 
the result that C vanishes in the low-energy regime ob¬ 
tained in the cases of random mass and random gauge 
potential will not be changed by higher order correc¬ 
tions. Nevertheless, the same conclusion cannot be sim¬ 
ply reached in the case of random chemical potential. 

In the spirit of RG theory [d8l-l50l|. the marginal na- 
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FIG. 1: Propagators: (a) fermion, (b) nematic, and (c) disor¬ 
der. 


ture of certain parameter is hardly stable against higher 
order corrections. It often happens that a marginal pa¬ 
rameter is turned to marginally relevant or marginally 
irrelevant once higher order corrections are taken into 
account. These two fates of disorder parameter may in¬ 
duce different behaviors of physical quantities. In or¬ 
der to address this issue, we need to go beyond the one- 
loop RG analysis, and make detailed two-loop (or even 
higher order) calculations. Interestingly, a recent work 
^ has revealed that two-loop corrections are able to 
make an important contribution to the quantum critical¬ 
ity of some particular three-dimensional Dirac semimet¬ 
als, which partly motivated the present work. 

In this paper, we combine the 1/fV-expansion method 
and the replica technique to analyze the interplay of ne¬ 
matic fluctuation and disorder scattering (in the follow¬ 
ing text we use ’’disorder” to uniquely denote random 
chemical potential). The perturbation expansion is car¬ 
ried out up to two-loop level. We derive a number of RG 
flow equations for fermion velocities vp and ua, velocity 
ratio va/vf, and disorder parameter These equations 
are coupled to each other, and thus need to be solved self- 
consistently. We show that the parameter C for random 
chemical potential becomes marginally relevant due to 
two-loop corrections. Therefore, the interaction between 
nodal fermions and random chemical potential is always 
in the strong coupling regime at nematic quantum crit¬ 
ical point, since any small parameter C flows eventually 
to infinity at the lowest energy. 

The rest of the paper is organized as follows. The effec¬ 
tive field theory and the corresponding Feynman rules are 
given in Sec. [TTl We present the RG transformations and 
perform detailed computations of one-loop correction to 
fermion self-energy and vertex functions in Sec. Ell which 
is followed by two-loop calculations in Sec.HV] In Sec.lVl 
we solve the RG equations and discuss the physical im¬ 
plications in Sec. IVII Finally, we briefly summarize our 
results in Sec. Em 


II. EFFECTIVE FIELD THEORY 


with being the standard Pauli matrices. The 

linear dispersion of nodal fermions originates from the 
(ia; 2 _j, 2 -wave symmetry of the superconducting gap of 
cuprates. Here, spinor tjj\ represents the fermionic quasi¬ 
particles (QPs) excited from nodal points (fif) and 
and the other two nodal points (isl - f^ . 
The repeated spin index a is summed from 1 to Nf, which 
is the number of fermion spin components. The ratio be¬ 
tween Fermi velocity vp and gap velocity va is rou^ly 
va/vf « 1/10, which is determined by experiments [l|. 

The action describes the Ising type nematic order 
parameter, which is expanded in real space as 



where r is imaginary time and c is velocity for cj). The 
mass parameter r tunes the nematic phase transition, and 
r = 0 defines the corresponding quantum critical point. 
Parameter uq is the strength of quartic self-interaction. 
The nematic order parameter couples to nodal fermions 
through the simple Yukawa term 


5^,0 = J d‘^xdT{Xo(l){tplaT^1pla + V’L'r"'V'2a)}- (4) 


The strong interaction between nodal fermions and 
nematic order can be handled by performing 1/iV^- 
expansion. The inverse of the free propagator of (j) be¬ 
haves as q'^+r. After taking into account the polarization 
function, an additional linear q-tevia will be generated, 
namely n(g) oc q. At the low-energy regime, the q-term 
dominates over the g^-term, which then can be neglected. 
Near the quantum critical point, we keep only the mass 
term and make the replacement that and 

r —>■ NfrXg, leading to 

S = (5) 

After integrating out fermion degrees of freedom, the ef¬ 
fective action for (j) becomes 


The interaction between massless nodal fermions and 
nematic order parameter is described by the following 
action [^. 

S = (1) 

where 5^ is the free action for nodal fermions 


The lowest-order Feynman diagram for the polarization 
function is shown in Fig. [2] and symbolizes the integral 

n(q,e) =j -^^^TT[T^Gl{k,u})T^G°^{k + q,uj + e)], 

where the free fermion propagator is 


S^ = 


d^k duj 
( 27 r )2 2tt 


-I- VpkxT'^ + VAkyT'^)i’la 


G°(k,a;) = 


ibJ -I- Vpk^T^ + VAkyT^ 


( 7 ) 
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FIG. 2: The polarization function for nematic order parame¬ 
ter. 


(a) 



FIG. 3: One-loop corrections to the fermion self-energy in the 
replica limit i? —>■ 0. 


As shown previously , the propagator for the nematic 
order parameter is given by 


G/{q,e) = n(q,e) 


1 (e^ + vUl) 


IGvfVa (e^ + Vpql + 

I 1 _ 

IGvfVa (e^ + Vpq^ + v'^q^yG 

in the vicinity of nematic quantum critical point r = 0. 
The Yukawa interaction has been studied extensively in 
recent years [l^. Huh and Sachdev found a stable fixed 
point at which the velocity ratio va/vf vanishes at the 
lowest energy. Such an extreme anisotropy then gives 
rise to a series of striking consequences, including anoma¬ 
lously scaling behaviors of specific heat, local density of 
states, and NMR relaxation rate 1/TiT [2^, reduction 
of thermal conductivity [2™, additional Cooper pairing 
between nodal fermions (28j |. strong suppression of su¬ 
perfluid density [27| and superconducting critical tem¬ 


perature and anomalous scaling of the penetration depth 


Now we introduce disorders into the above model. Dis¬ 
orders exist in almost all realistic condensed matter sys¬ 
tems and are known to be responsible for many of the 
low temperature properties. In the present problem, the 
nodal fermions may interact with three sorts of disorders, 
which can be represented by the following general term 


5'dis = J (fx.yyx.)Ty{x)A{x.), 


(9) 


where the matrix F is F = I for random chemical poten¬ 
tial. It will be replaced by in the case of random mass 
and random gauge field. It is traditional to assume 
that A(x) is a quenched, Gaussian white noise potential 
defined by the following correlation functions 

(A(x))=0; (A(xi)A(x 2 )) = C<5^(xi - X 2 ). (10) 

The disordered potential A(x) is randomly distributed 
in space, and needs to be properly averaged. A commonly 
utilized method to av erag e disorders is to introduce the 
so-called replica trick j^, , which states that one 

can replicate the partition function Z hy R times with R 
being a positive integer, Z —>• Z^. After doing so, one 
obtains a useful identity 

Z^-1 

\nZ = lim -. (11) 

R^o R ^ ’ 


Though R is initially supposed to be an integer, we can 
regard it as a continuous variable and simply take the 
limit i? —>■ 0. It is now technically feasible to average 
over the replicated partition as follows 

WP[A]Z^, (12) 

where 

P[A] = exp (^-J y d"xA2(x)^ . (13) 

We then insert the partition function Z and integrating 
out the random potential A(x), and finally obtain an 
effective action 

^ + + d^drdr' 

XV’i(x, T)?/'m(x, t)-!/'1(x, r')V'„(x, t'), (14) 



where m, n are the replica indices. To proceed, it is more 
convenient to rewrite the disorder-related term in the mo¬ 
mentum space. Therefore, we have 


SeS 


nm I nm , nm 


X 


C j" d^lid^h'd^h"dujdu! 

lj W 

V’i(k",a;")V'n(k + k" - 


-^m(k,w)'0m(k',w) 

k',c.")- (15) 


The propagators of the replicated action are delineated 
in the Fig. [T] There are two noticeable features of this 
effective action: (i) the scattering of fermions by static 
disorder does not exchange energy and hence the energy 
of nodal fermions is conserved; (ii) the quartic interaction 
connects different replica indices. 

We here list the one-loop and two-loop Feynman di¬ 
agrams contributing to the fermion velocities and dis¬ 
order strength: (i) There are two one-loop diagrams in 
the replica limit i? —>■ 0 which contribute to the fermion 
self-energy as represented in Fig. [31 (ii) There are four 
one-loop diagrams in the replica limit i? —>■ 0 which con¬ 
tribute to the disorder vertex as represented in Fig. jj] 
(iii) There are seven two-loop diagrams in the replica 
limit i? —>■ 0 which contribute to the fermion self-energy 
as represented in Fig. |S] (iv) There are thirty two-loop 
diagrams in the replica limit i? —>■ 0 which contribute 
to the disorder vertex as represented in Fig. |6l In the 
next section, we will calculate all the associated Feyn¬ 
man diagrams and derive the RG flow equations for all 
the physical parameters appearing in the effective action. 
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(a) (b) 





(c) 


(d) 


FIG. 4: One-loop corrections to the disorder vertex in the 
replica limit i? —>■ 0. 


III. ONE-LOOP RG ANALYSIS 

In this section, we first make RG analysis at the one- 
loop level. Different from Ref. [i^, the disorders are 
averaged here by using the replica method. The general 
scheme presented in this section will be applied directly 
to perform two-loop calculations in the next section. 

In order to do RG calculations, we make the following 


scaling transformations [l^, IH, , 

h = k[b, (16) 

oj = oj'b, (17) 

^i.2(k,a;) = (18) 

<^(q,e)=0'(q',e')e5/o(5-'7^)'^', (19) 


where I is a freely running length scale with b = e~^ and 
i = x,y. The scaling parameters rif and rib will be de¬ 
termined by the self-energy and fermion-nematic vertex 
corrections. Notice that the energy is required to scale in 
the same way as momentum, which means the fermion 
velocities are forced to flow under RG transformations. 

With the help of Dyson equation, we know that inter¬ 
actions induce a self-energy correction to the free fermion 
propagator, i.e., 

G^^(k,a;) = -iuj+vpkxT''+VAkyT'" (20) 

where E(k, w) is the fermion self-energy function. In 
the current problem, E(k, w) receives corrections from 
both fermion-nematic interaction and fermion-disorder 
interaction. At the one-loop level, formally we have 
E(k,u;) = .sphere and 

j^Li(b) generated by nematic order and disorder, re¬ 
spectively. 


A. Fermion self-energy 

The one-loop diagrams for fermion self-energy are de¬ 
picted in Fig. [21 The one-loop nematic self-energy, shown 


in Fig. jSKa), has already been calculated earlier by Huh 
and Sachdev who found that 


jyLl(a)/'n 

- = + C2VFk,T^ + C3VAkyT\{21) 

am A 

where the coefficients C\, C 2 , and C 3 can be found in 
theO The one-loop disorder-induced fermion self-energy 
shown in Fig.jSKb), is given by [53| 







rG°(k,ca)r 


() factor 
4 2'kvfva 


iui In A. 


( 22 ) 


This expression tells us that is independent 

of the concrete form of the matrix F. Another impor¬ 
tant feature is that E^^*^*’)(ia;) is independent of momen¬ 
tum, which reflects the fact that the quenched disorder 
is static. It is now easy to get 


dELi(b)(^^) 
din A 


= 2 Cgiuj, 


where 


Ca = 


C 1 

4 2fvfva 


(23) 


(24) 


B. Disorder vertex 


In the replica limit i? —>■ 0, the one-loop diagrams for 
the corrections to fermion-disorder vertex are given in 
Fig. m The diagram of Fig. |4Ka) can be calculated by 
emp loying the method proposed by Huh and Sachdev 
|24l| . At zero external momenta and frequency, the cor¬ 
responding vertex correction is expressed as 



where Q = (q, e) is a three-momenta. Here, /C(y) is an 
arbitrary function with /C(0) = 1, and it falls off rapidly 
with y, e.g., /C(y) = [^. However, the results are 

independent of the particular choices of lC{y). The above 
equation can be converted to 


dH 


dis 


din A 


C\ 

4 / Stt^ 


dx / d9H{Q), (26) 


where 


1 


HiQ) = 


1 


Nf {—ivFX + vf cos9t^ + vasyuOt^) 
1 

{—ivpx + vf cos 9 t ^ + va sindr®) 

. 1 


I 


XT 


n(g) 


(27) 
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Here, the matrix I corresponds to the coupling between 
nodal fermions and random chemical potential. After 
straightforward computation, we have 


din A 



(28) 


where 

^5 


dx 


d0 


2(va/vf) 

Nfir^ j-oo JO 
(x^ — cos^ 0 — (va/vf)'^ sin^ 0) 
(x"^ + cos^ 0 + (va/vf)"^ sin^ 0)"^ 


G{x,0) 


= -Cf 


(29) 


Fig. IlKb) is the vertex correction due to disorder averag¬ 
ing, and given by 



/ i ‘2 

^IG°(cc,p)IGO(cc,p + k)I. (30) 


Taking external momentum k = 0 and keeping only the 
leading divergent term, we find that 


where 




Ll(b) 


dis 


din A 



Gr 


C l c. 

4 2fvfva ^ 


(31) 


(32) 


The contributions from Fig. 0Kc) and Fig. HKd) cancel 
each other M- 


C. RG Equations with disorder at one-loop level 

RG theory requires that the momentum shell between 
6 A and A should be integrated out, while keeping the 


—iuj term invariant. Using the one-loop contribution to 
fermion self-energy, we have 


rbA 


d^kdcij'01 


—iu! — Ci(—iui) In — 


+ 2 Cg{-iuj)\n^ 


iP 


/ bA. 

d’^kdujip^(-iw) [1 -h ( 2 Gg - Ci)l]i; 






(33) 


After scaling transformation, this term goes back to the 
corresponding free form, which implies that 

Vf=2Cg-Ci, (34) 

at the one-loop level. The kinetic terms should also be 
invariant under scaling transformation, so we obtain the 
following RG equations up to one-loop level for fermion 
velocities: 


dvF 

~dr 

dVA 

~dr 


(Gi -C 2 - 2Cg)VF, 

(35) 

(Cl - G 3 - 2 Gg)^;A. 

(36) 


Based on these expressions, the one-loop RG equation for 
the ratio between gap velocity and Fermi velocity can be 
derived as 


-^ = (C2-C3) — . (37) 

dl vf 

The disorder parameter ^ indirectly appears in the 
above equations which exists in the Cg. Due to the in¬ 
terplay of nematic fluctuation and disorder scattering, 
this parameter also flows under RG transformation. In¬ 
cluding one-loop correction due to nematic and disorder 
interactions yields 


J d‘^kd‘^k'd‘^k''dujdLo''^p)n(k,u})'iprn(k',u}'')'ipn(k + 'k'' - k '~ ^ 

J d2kd2k'd2k"da;dw'Vi(k,w)V'™(k',w)V't(k",w")V'n(k + k"-k',w") (^-0 I[l + 2(2Gg-Gi);] 

J d^kd^k'd^k"du}du}"'ijjl^{k,uj)tprn(k',uj)tpl^(k'',uj")tpn{k + k" — k',uj") ^0 _ ( 33 ) 


After redefining energy, momentum, and field operators, we have 

rA 


J d^kd^k'd^k''du}du}"ipj^(k,u;)'iprn{k',uj)ipn{k",u;'')ipnik+ k" — k',ijj") (39) 
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Since rjf = 2Cg — Ci, it is easy to obtain the following 
one-loop RG equation for 


From this result, we know that the strength parameter 
for random chemical potential C is marginal at the one- 
loop level. 

IV. TWO-LOOP RG ANALYSIS 

As shown previously in Ref. [i^ and also in the last 
section, we have found that the strength parameter of 
random chemical potential is marginal at the one-loop 
level near the nematic quantum critical point located in 
a d-superconductor. According to the RG theory, the 
fate of one marginal parameter can be changed once 
higher order corrections are incorporated. Specifically, 
the marginal parameter may be turned to either rele¬ 
vant or irrelevant by the corrections. Recent work found 
that the two-loop corrections do make a significant con¬ 
tribution to the quantum criticality of three-dimensional 
Dirac semimetals M- It is therefore important to ex¬ 
amine the fate of marginal random chemical potential 
against higher order corrections. To address this issue, 
we now go beyond the one-loop analysis given in Ref. 
and perform a detailed two-loop RG analysis. 

The Feynman diagrams at two-loop level are presented 
in Fig. [5] and Fig.[S] These diagrams can by calculated by 
employing the methods introduced in [ 2 ^, [s^, and 
1^ . Paralleling the procedures used in Sec. IIIIl we will 


be allowed to derive the RG equations at two-loop level. 
All the related coefficients are demonstrated in the IS (i) 
mi, 1712, and m 3 collect the contribution from Fig. [5] (a); 

(ii) di, d 2 , and da collect the contribution from Fig.[S](c); 

(iii) /i, /a, and /a collect the contribution from Fig.[5](e) 
and by comparing with results with Ci, C 2 , and C 3 , we 
can obtain fi = Ci, {i = 1,2,3); (iv) Cu collects the 
contribution from Fig. |n](g2); (v) C/ collects the contri¬ 
bution from Fig.[6](g3); (vi) Ck collects the contribution 
from Fig. [5] (i2); (vii) C; collects the contribution from 
Fig-ini(j4); (viii) Cm collects the contribution from Fig. [5] 
(k2); (ix) Cn collects the contribution from Fig. |n](k3); 
(x) Others Feynman diagrams in Figs. [5] and | 6 ] either 
counteract each other due to the appearance of ladder 
and crossing diagrams or their contributions can be ex¬ 
pressed by Cl, {i = 1,2,3), Cg and/or above two-loop 
coefficients. 

The Dyson equation for fermion propagator can now 
be formally written as 

G/;^(k,a;) = -ioj + VFk,j;T''+ VAkyr"^ 

-S^i(^)(k,a;)-E^i('")(k,w) 

-I]^2('=)(k,a;) - 

Before computing the RG equations, we need to obtain 
the expression oi rjf. Integrating over the momentum 
shell between bA and A, we obtain 


J d2kda;^/1'(-iw)|l-(Gi-2Gg)Z- 


mi 


C2 Ga 


-^-Cgdl+Cg^Cl-^-^]- 2Cgfl - 2C^g - 2Cg 


u 


J d‘^kduijj\-iu}) - {Cl - 2Cg)l- 


J d^\iduj'tp'^{—iu})expl{2Cg — Ci)l — 


^ - Ggdi + Gg (^Gi - ^ ^ - 2Gg/i - 4Cf 
^ - Ggdl + Gg fci - ^ - 2Gg/l - 4G/ 


AlP, 


(42) 


where mi and /i come from the contribution of Fig. [5] Therefore, rjf can be obtained at the two-loop level, 
(a) 


After scaling transformation, this term should go back 
to the free form, thus 


rjf = (2Gg - Gi) - + Gg ( 2 G 1 -C 2 -C 3 - 2di 


- 4/1 - 8Gg) 


(44) 


J lJfdl+\^{2Cg-Cl)l- 


mi 


-G„d 


ff«i 


+ Cg[Cl-^-^]-2Cgfl-ACf 


P ^ = 0. (43) 


A. RG equations of fermion velocities 

After renormalization and rescaling transformation, 
the Fermi velocity term vpkxT^ must return to its origi- 
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FIG. 7: The running of va (left) and v_f (right) without disorder at the two-loop level. 


X exp < — 


^-( 2 Cg-Cl)l + 


mi 


-C„di 


+Cg ^ - 2Cgh - 4C|] 

xexpj-CaZ- ^ - ^3(^2 + 2/2) 



d? {v'pkxT^)^^ ^ 


( 45 ) 


where the renormalized Fermi velocity v'p becomes scale- 
dependent and is determined by the following expression 


v'p = VF exp |(Ci - C2 - 2 Cg)l + g + Cg (Ci - 


(^2 

T 


-^-di-2fi + {d2 + 2f2)-4Cg 


( 46 ) 


It is now easy to extract the RG equation for vp up to two-loop level, 


^ ^ (Cl - C 2 - 2 C 3 ) + 


3 + 2 Cg (2C2 - Cl + d2 - di - 4C3) - C3 (C2 + C3) 


Ij vp, ( 47 ) 

where the identity fi = Ci is used. The two-loop RG equation for the gap velocity va can be analogously obtained, 
i.e., 

^ = |(Ci-C3-2C3)+ ~ + 2 Cg (2C3 - Cl + d3 - rfi - 4C3) - C3(C2 + C3) l^VA- ( 48 ) 

Based on the above two equations, we can also derive the RG equation for the velocity ratio up to two-loop level, 

~ + 2 Cg (2C3 - 2C2 + da - d2)] 4 — ■ ( 49 ) 

o \ ) '<^F 


d^ 


B. RG equation of disorder strength parameter 


The strength parameter vp of random chemical poten¬ 
tial enters into the above equations for fermion velocities. 
Due to the interplay of nematic fluctuation and disorder 
scattering, this parameter also flows under RG transfor¬ 


r 


mation. We calculate the two-loop diagrams for the cor¬ 
rection to fermion-disorder vertex function and combine 
them with one-loop contribution, and then obtain 


rbA 


d^hd^h'd'^k"du!du}"'il)l^(\i, w)'^/'m(k^ w)'!/'!(k", w'‘ 
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l l 

FIG. 9: Flows of the disorder strength with the initial values ( va / vf )^ = 0.1, = 0.1, = 10“® (left) and 

= 10”“^ (right) at the two-loop level. 



j cP'kcfk'd^k"dwduj" w)^m(k', uj)'tpl^{k'', uj") 

xip„{k + k'' - k', w") ^-0 exp| - 2Cgdi 
+12C^g + CiCh + CgCf + 2Cf - SCiCg + Ck + Ci 
+Cm + CgCr, + Cg(C2 + Cs)] /'} 

pK 

/ d^kc?^k'(i^k”da;(ia;"'0)!)j(k, a;)'!/'m(k^ UJ)'^jJ'l^{k", oj") 


xipnikd- k" - k',uj") ( - 


c:' 

4 ; ’ 


( 50 ) 


where 


FIG. 8: The running of va/vf without disorder at the two- 
loop level. 


x-ipnO^ + k" - k',a;") - 


c 


1 - 2Cil + ACgl + ACn 


2,2 


f 

-f2,2 


g 

2,2 


+Cr + CiChl" + CgCfl^ + Ckl^ - 2CgCiV + 


+2Cil^ + Cml" - ACiCgV - Sqr + + CgCnr 

+2Cg{C2+Cs)l^'^ 



2Cgdl + 12 C| + ClCh + CgCf 


+2Cl - SClCg + Ck + Ci + Cm + CgCn 


+ Cg{C2 + Cg) 



( 51 ) 


The RG equation for parameter C, up to two-loop level is 
finally found to be 


1 


dl 


21 


CgiUCg - 2dl + Cf) + Cl(C,, + 2 Cl - SCg) + 

I— 


( rni 

V 3 


+ Ck+Cl+ Cm) + CgiCr, +C2+ Cg)] C- ( 52 ) 


V. NUMERICAL SOLUTIONS FOR RG 
EQUATIONS AND DISCUSSIONS 

The RG equations for fermion velocities and disorder 
parameter have been derived in the last section. The aim 


of this section is to analyze the solutions of these equa¬ 
tions, which are self-consistently coupled to each other 
and therefore need to be solved numerically. To make the 
discussion more transparent, we first consider the clean 
limit with C = 0 and then include disorder later. 
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FIG. 10: Flows of ^ with initial values {va/vf)° = 0.1, = 1, = 0.1, = 10“® (left) and = 10“^ (right) at the two-loop 

level. 


A. Clean limit 


In the clean case, C = 0, the RG equations have the 
following form up to two-loop level, 


dvp 

~dr 

dVA 

~dr 


di 


C1-C2 + 
Cl - C3 + 
C2 - C3 + 


(mi - m2) 


(mi - m 3 ) 


(m 2 - m 3 ) 


VF, 

VA, 

VA 

VF 


(53) 

(54) 

(55) 


Solving these RG equations numerically, we find that 
the results qualitatively agrees with those obtained at 
one-loop level Ref. [l^, which can be clearly seen from 
Figs. © and ([5]). In particular, vp vanishes with grow¬ 
ing I much more slowly than va, so the velocity ratio 
va/vf ^ 0 ai large length scale, which corresponding to 
a stable fixed point of extreme anisotropy. It therefore 
turns out that the one-loop results are robust and that 
the 1/iV-expansion is fairly reliable. 


B. In the presence of disorders 

Now we consider the influence of random chemical po¬ 
tential. Before running into complex numerical compu¬ 
tation, it is helpful to first make a simple analysis about 
the behavior of parameter (). From the flow equation 
of c dSH), it is convenient to extract a general differen¬ 
tial equation, ^ = CIC- Apparently, the fate of C in 
the low-energy region (i.e., large I limit) is unambigu¬ 
ously determined by the value of the coefficient C: as 
the running length scale I increases, C, flows towards in¬ 
finity, vanishes, and stays at a certain constant for C > 0, 
C < 0, and (7 = 0, respectively. Previous one-loop cal¬ 
culations [ 4 ^ showed that the coefficient C vanishes, 
so the strength parameter C of random chemical poten¬ 
tial is marginal and remains a constant as I varies. After 


including two-loop corrections, C may become either pos¬ 
itive or negative, which will be addressed by solving the 
coupled flow equations gzn, (gg), dUl), and (El]). 

The numerical solutions of RG equations with prop¬ 
erly chosen initial values of disorder parameter are shown 
in Figs, m and [TOj As I grows, C increases quickly and 
eventually diverges in the limit I —>■ -l-oo. This behavior 
does not depend on the initial values of the running pa¬ 
rameters. Therefore, the strength parameter of random 
chemical potential, which is marginal at one-loop level, 
becomes relevant after including the two-loop order cor¬ 
rections. This change may have significant influence on 
the physical properties of the system under consideration 
since even an infinitesimal relevant parameter will flow to 
very large at ultra low energies. 

A natural question arises: how can we understand the 
marginally relevant nature of the characteristic parame¬ 
ter C of random chemical potential and its impacts on the 
nematic quantum criticality in d-wave superconductors? 
In the clean limit, C = 0, the strong quantum fluctuation 
of nematic order results in non-Fermi liquid behavior of 
nodal fermions [13, HO] . Once random chemical potential 
is introduced, the non-Fermi liquid state is immediately 
turned into a diffusive metallic state which has 

a finite zero-energy density of states and finite scattering 
rate. This conclusion is valid for all possible values of 
since irrespective of its initial value, unavoidably flows 
to infinity in the low-energy regime. 


VI. PHYSICAL IMPLICATIONS 

In addition to the tendency of non-Fermi liquid behav¬ 
ior due to the marginally relevant behavior of random 
chemical potential, we, within this section, try to de¬ 
lineate investigate the influence of nematic fluctuation, 
gapless QPs and marginally relevant random chemical 
potential on a number of significantly physical observ¬ 
ables in a d-wave cuprate superconductor at the nematic 
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QCP, such as the superfluid density, critical temperature, 
and thermal conductivity. 

Learning from the Fig. |9]and the discussion in Sec. El 
we can be recalled some points. On the one hand, 
the strength of random chemical potential, due to the 
marginally relevant behavior, flows away with the energy 
scale I reaching a critical value, dubbed I*, at the nematic 
QCP. This indirectly triggers the unphyscial behaviors of 
fermion velocities vp and va at a very scale 1 = 1* which 
is dependent on the initial values of vp and va by the 
growing strength of random chemical potential and will 
be carefully investigated in the following subsections. On 
the other hand, runaway behavior of physical quantities 
suggests a first-order phase transition [bTH^ . In order 
to figuratively describe it in finite t emp eratures, we em¬ 
ploy a useful formula T = TcC”* (Tc is the critical 

temperature of superconductor) to translate I* to T*. 

It is necessary to further emphasize that our effec¬ 
tive theory and method are only valid for the continu¬ 
ous phase transitions and hence we subsequently focus 
on the critical behaviors of physical observables with the 
confined region T > T*, such as superfluid density, criti¬ 
cal temperature and thermal conductivity. 



FIG. 11: Temperature-dependence superfluid deusity iu the 
preseuce of the nematic fluctuation and random chemical po¬ 
tential at nematic QCP for the representatively initial values 
of va/vf and random chemical potential. The primary con¬ 
clusions are independent of the initial values. 


A. Superfluid density and critical temperature 


As studied by Lee and Wen , the superfluid density 
in the noninteracting case exhibits a linear temperature 
dependence which is in agreement with experiments 741 


Ps{T) _ PsjO) _ 2\n2vp ^ 
m m TT VA ' 


(56) 


with Ps(0) = ^ where x and a represent the dop¬ 
ing concentration and lattice spacing, respectively and 
Tc (X by defining Ps{Tc) = 0, reproducing the 

Uemura plot [t^. The other term of right hand side 
is the contribution of the normal QPs density. As pre¬ 
sented in the previous section, the fermion velocities vp, 
va, and the ratio va/vf are heavily renormalized by the 
interplay between nematic fluctuation and the effects of 
random chemical potential in the vicinity of the nematic 
QCP. We can approximately derive the superfluid den¬ 
sity of of the d-wave superconducting state at the nematic 
QCP which is renormalized by nematic fluctuation and 
marginally relevant random chemical potential , 

p',{T)=pM-pniT), (57) 


are dependent on the running of velocities vp and va 
and also the behavior of disorder which is coupled with 
the velocities vp and va as lineated by Eqs. dUl), dig), 
and (H51) . To include these corrections and after some 
numerical calculations, we obtain the p(,(T)/ps(T) with 
the different initial values of fermion velocities vp and va 
for T > T* as depicted in Fig. [TT] 

Studying from the numerical results in Fig. 1111 we can 
draw a conclusion that the superfluid density and crit¬ 
ical temperature are largely suppressed in the presence 
of marginally relevant random chemical potential nearby 
the nematic QCP. This may be very instructive to pro¬ 
vide another clue to locate the nematic QCP. 


B. Thermal conductivity 


We then turn to the thermal conductivity of QPS. By 
paralleling the derivations in Ref. [t^, we obtain the 
thermal conductivity of QPs for the constant values of 
fermion velocities vp and ua, 


K 



T ~ 

3 

<1 


(59) 


where the normal QPs density is complicatedly depen¬ 
dent on the fermion velocities vp and va 

Pn{T) _ 4 r Vp{k)e >‘bt 

m kpT J (27r)2 ^ 

\l + e >‘bt 



and for the energy scale-dependence vp and va which are 
comprised the nematic fluctuation and marginally rele¬ 
vant random chemical potential in the vicinitv of nematic 
QCP, 

T " [tjJo (TTW 


U_F(k) UA(k) 


UA(k) Vp(k) 


,(60) 


By combing the Eqs. (EZD and (1551) and Eqs. (H71) . (H51) . with Ag determined by the lattice constant and Co the 
and (gg, the superfluid density and critical temperature initial value of disorder strength. In order to examine 
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FIG. 12: Thermal conductivity of QPS in the presence of ne¬ 
matic fluctuation and random chemical potential at nematic 
QCP for the representatively initial values of va/vf and ran¬ 
dom chemical potential. The primary conclusions are inde¬ 
pendent of the initial values. 


the influence of nematic fluctuation and random chemical 
potential, we numerically carry out the Eqs. (03), 03, 
03, 03 and 03 and reach the numerical results as 
presented for T > T* in Fig. [T3 

As exhibited in Fig. the thermal conductivity of 
QPs, influenced by nematic fluctuation and marginally 
relevant random chemical potential at the nematic QCP, 
is highly increased and insensitive to the initial values 
of random chemical potential. This may be helpful to 
enhance our understand the phase diagram of d-wave su¬ 
perconductors. 

VII. SUMMARY 

In summary, we have investigated the impact of a par¬ 
ticular disorder, namely random chemical potential, on 

I 


the stability of nematic quantum critical point located 
in the superconducting dome of d-wave superconductors. 
We performed a detailed RG analysis up to two-loop or¬ 
der within an effective field theory that describes the in¬ 
terplay of nematic fluctuation and disorder scattering. 
The effective parameter that characterize the strength of 
random chemical potential is marginal at one-loop order, 
but becomes relevant due to the inclusion of two-loop or¬ 
der corrections. This means that, even if we start from 
a very weak disorder, the effective strength of disorder 
eventually becomes infinitely large at the lowest energy. 
The Dirac fermions therefore enter into a diffusive metal¬ 
lic state driven by random chemical potential at the ne¬ 
matic quantum critical point [571 l6lM65|| . Furthermore, 
we carefully study the critical behaviors for a number of 
significantly physical observables in a d-wave cuprate su¬ 
perconductor at the nematic QCP, such as the superfluid 
density, critical temperature, and thermal conductivity 
in the vicinity of nematic QCP. 
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Appendix A: Some coefficients used in the main text 


The coefficients for one-loop Feynman diagrams are 


Cl 


C2 


C 3 


2va 


NfTT^Vp 


^2ii- — cos^ 9 — 

dx / dO 




sin^ 0 


-00 Jo 


x"^ + cos^ 9 -\- 






sin"^ 9 


NfTT-ivF d_oo Jo 


■9-x^- (^^^in^d 


x'^ + cos^ 9 -I- 


2va r. ^ “ 

dx / dO - 




sin^ 9 


■G(x,9), 


sin^ 9 


NfTT^Vp 


x"^ + cos^ 9 -I- sin^ 9 


-^g{x,9), 


(Al) 


(A2) 


— 00 


0 


(A3) 
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= 


+ cos^ ( 


+ sin^ ( 


+ sin^ 0 Jx'^ + sin^ 9 + cos^ 9 


C ^ ^ . 

® 4 27rt!i7’fA 

For two-loop Feynman diagrams, the corresponding coefficients are 


^oo /*27r 


Gix,e) 


x^ -|- cos^ 6 + 


mi = — —r— / dx d 6 - 

Nfn^F J_^ Jo 


xUx^ Cix^ + { 2 C 2 - Cl) cos"^ e + {2C3 - Cl) {^j sin^e 


- iCix"^ + (2C2 - Cl) cos^ 9 + {2C3 - Cl) f ^ j sin^ 9 


X x"^ + cos^ 0 + ( — I sin^ 9 


rco /*27r 


m 2 = ——^— / dx d9' 
Nf^Hp y_oo io 


G{x,9) 


x'^ + cos^ 9 + sin^ 9 


x| {C 2 - 2 Ci)x^+ {C 2 - 2 C 3 )(^j sin^ 6» - 3(72 cos^ 0 
/ \2 1 

X -I- cos^ 0 -I- ( — ) sin^ 0 — 4cos^0 ((72 — 2Ci)x^ 

+{C2 - 2C3) (—^ sin^ 9-C2 cos^ 9] \ , 

\VF J J J 


nOO /*27r 


m 3 = — — 3 — / dx d9' 
NfTT'^VF J-oc Jo 




Gix,9) 


x'^ + cos^ 9 + sin^ 9 


X x^ -I- cos^ 9 + (—^ sin^ 9 — 4 (—^ sin^ 


sin^ 9 ( 2 Ci — C3)x^ 


ffiCsf—1 sin^^+(2(72-t(73)cos2 6ij 1 
\VF/ J J 


^ r ,, Cm _ cci _^ 

° x2+cos2 0+(^)%in2 0 

x/ 3x^ — cos^ 0 — ^ sin^ 0 x^-I-cos^ 0-I- ^ ^ sin0^ 


-4x2 x2-cos20- (^ ) Sin20 , 
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with 



(AlO) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 
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E2 


Va 

+ cos^ ^ — 1 — I sin^ (p 

Vp 


Va 

+ 1 — I sin^ p + cos^ p 

Vp 


-2j/" 


+ cos^ iy9 — 3 


VF 


sin^ ip 


^y/2v^ f°° {cos^ip-y‘^)Q{y,ip) 

iw.L ‘‘’'I - 


p J —oo Jo 


8V2'. 


VA 


oOO n2'K 


J_ 


dy / dip- 


2/2 + cos2 if + sin"^ if 

Qiy,^) 


n3’ 


p J —oo Jo 


VA 


— I sin^ ifi-y^ 
VF 


2/2 + cos2 ip + sin^ if 


Va 

— I sin^ ip + cos^ tp + 2/^ 
VF, 


+2 


VA 

VF 


■ 2 

Sin ip 


3y^ — cos^ ip — 


VA 

VF 


■ 2 

Sin ip 


(A18) 

(A19) 


(A20) 
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